Clustering properties of a generalised critical Euclidean network 
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Many real-world networks exhibit scale-free feature, have a small diameter and a high clustering 
tendency. We have studied the properties of a growing network, which has all these features, in 
which an incoming node is connected to its ith predecessor of degree ki with a link of length £ 
using a probability proportional to kfl a . For a > —0.5, the network is scale free at B = 1 with 
the degree distribution P(k) oc k ' and 7 = 3.0 as in the Barabasi- Albert model (a = 0, 8 = 1). 
We find a phase boundary in the a — 8 plane along which the network is scale-free. Interestingly, 
we find scale-free behaviour even for 8 > 1 for a < —0.5 where the existence of a new universality 
class is indicated from the behaviour of the degree distribution and the clustering coefficients. The 
network has a small diameter in the entire scale- free region. The clustering coefficients emulate the 
behaviour of most real networks for increasing negative values of a on the phase boundary. 
PACS numbers: 05.70.Jk, 64.60.Fr, 74.40. Cx 



Recent studies of many complex real-world networks 
of diverse nature, e.g., social networks, biological net- 
works, electronic communication networks etc. reveal 
some striking similarities in their underlying structures 
[1]. The diameter T>, a measure of the topological exten- 
sion of the network, the clustering coefficient C, a measure 
of the local correlations among the links of the network 
and the nodal degree distribution P(k) are some of the 
few important quantities which exhibit the similarities 
among the different networks. Many of these networks 
exhibit small- world network (SWN)-like properties [2], 
i.e., the diameter V(N) of the network scales logarithmi- 
cally with the number of nodes N while the clustering 
coefficient has a high value. In some of these networks, 
there is no characteristic scale manifested by the typi- 
cal power law decay of the tail of the degree distribution: 
P(k) oc & -7 [3], where P(k) is the number of nodes which 
are linked with k other nodes. These networks are called 
scale- free networks (SFN). 

Typically, the clustering coefficient measures the con- 
ditional probability that an arbitrary pair of nodes are 
linked, provided both are linked to a third node. The 
clustering coefficient can be studied as a function of two 
different variables: C(N), the clustering coefficient per 
node averaged over all N nodes as a function of the net- 
work size N and C(k), the clustering coefficient per node 
averaged over all nodes with degree k as a function of k. 
Obviously C(N) = Y, k P(k)C{k)/Z k P{k). 

In some recent studies [4,5], it was shown that several 
real networks, like the actor network, language network, 
the Internet at the autonomous system level etc., which 
are known to exhibit scale-free behaviour and have small 
diameters, have another common feature, i.e., C(k) has 
a power law dependence: C{k) oc k~ x whereas the total 
clustering cocffcient C(N) has a high value. 

Attempts to capture the three features of small di- 
ameter, high clustering and absence of a characteristic 
scale, which occur in many real world networks, in a sin- 
gle model, have been faced with certain difficulties. The 



first model to mimic a small- world network is the Watts- 
Strogatz model (WS) [2]. Here the nodes are arranged 
on a ring with links to the nearest neighbours and small- 
world features can be achieved by re-wiring the nearest 
neighbour bonds to randomly link an arbitrary pair of 
nodes even with a very small probability. However the 
nodal degree distribution in the WS model failed to show 
scale-free feature. The Barabasi-Albert (BA) model is a 
prototype for a SFN in which the network is grown by 
adding nodes one by one, and a new node gets attached 
to an older one with a probability proportional to its de- 
gree. Although the scale-free property was successfully 
achieved and the network had a small diameter, the clus- 
tering coefficient C(N) showed a power law decay with N 
(C(N) oc N~ '' 5 ), while C(k) remained a constant with k 
[1,4], thus failing to capture the feature of high clustering 
tendency of real networks. 

Successful attempts to capture all the desirable fea- 
tures of a network have been made by defining other 
models [4,6-10] subsequently. For example, in a deter- 
ministic growing graph [6], which is argued to simulate 
a citation network, exact calculations showed that it has 
small diameter, scale-free feature as well as C(k) oc 1/k. 
In [8,10], suitable modifications are done to generate tri- 
ads (and consequently a high clustering coeffcient) in an 
otherwise BA type of growing network. In [7], an old 
node is deactivated with a probability proportional to 
its inverse degree in a growing network to get a high 
clustering coefficient. In [9], spatial distances have been 
incorporated in some specified manner which also gave 
the desired features of a real network to a large extent. 
A power law dependence of C(k) can be obtained in de- 
terministic and stochastic scale-free networks with hier- 
archical structure also [4]. 

While in a majority of real- world scale- free networks 
C(k) oc 1/k , some other networks like the Internet router 
network, the power grid network [4,5] of the Western 
United states and the Indian railway network [11] (which 
does not have scale-free behaviour) showed a different 
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behaviour: C{k) shows no dependence or logarithmic de- 
pendence on k. In [4], this behaviour was argued to be 
due to the presence of geographical organisation in such 
networks in the sense that there are actual physical con- 
nection between the nodes and the networks are defined 
in real space. A comparison of the clustering coefficients 
in a model network with and without geographical organ- 
isation could therefore help to understand the relation 
between geographical organisation and clustering better. 
It should be pointed out here that C(k) is also a constant 
in the BA model where a metric is not defined. 

In a network defined in real space, the spatial distance 
between the nodes is expected to play an important role 
in constructing the links. On the other hand, the rule 
of preferential attachment has been very successful in 
achieving the scale-free feature and small diameter of a 
network. We have therefore considered a growing net- 
work in which both the preferential attachment as well 
as the spatial distances are parametrically incorporated 
in the attachment probability and can be independently 
tuned. Here we would like to mention that although some 
networks are not defined in real space, spatial distances 
are still expected to be implicitly involved, e.g., in a so- 
cial network, people in the same locality are much more 
likely to know and influence each other. Although the 
concept of geographical locality does not exist explicitly 
in all networks, one can still define a "closeness" factor 
in many networks, e.g., in the citation network, a paper 
is likely to be cited with a higher probability when the 
contents of it is "close" to that of the citing paper. 

The network we have under consideration is evolved 
from the time t = and at time t an incoming node gets 
attached to the i-th node with degree hi, at a distance £, 
according to the following probability: 

^(t) ~ fcf (t)r. (i) 
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FIG. 1. The phase diagram of the network in the a — /3 
plane. Scale-free behaviour is observed only at the boundary. 
The point A indicates a change in the critical behaviour: to 
the right of A the critical behaviour is of the BA type while 
to its left we find a new critical behaviour. 



In some earlier studies [9,12-15], spatial dependence 
in a growing network has been studied where the attach- 
ment probability is dependent on the spatial distances 
between the node. In [9], the clustering coefficients were 
also calculated. However, the spatial dependence was not 
incorporated in a way it could be systemetically studied. 
In the present study, the aims are (a) to identify the re- 
gions where the network is scale free in the a — [3 plane 
and (b) to study the behaviour of the clustering coeffi- 
cient as a function of the parameters a and (3 and (c) 
to check whether the diameter of the network and the 
average shortest distances scale logarithmically with the 
number of nodes. 

In this network, each incoming node gets bonded to 
to distinct nodes. In order to study clustering properties 
to should be at least equal to two (m = 1 would lead 
to a tree like structure with no loops and all clustering 
coefficients are trivially zero here.) Results for some lim- 
iting cases of the model defined by (1) are known. The 
a = and (3=1 case corresponds to the scale-free BA 
network [3]. Networks with a — and arbitrary values 
of (3, considered in [16], showed that scale- free behaviour 
existed only for (3 = 1. For (3 > 1, there is a tendency 
of the incoming nodes to get connected to a single node 
and this behaviour is termed "gelation". For (3 < 1.0, 
the behaviour of the degree distribution is stretched ex- 
ponential. The effect of Euclidean distances were incor- 
porated in a BA kind of network [13,14] by keeping a 
non-zero and (3 = 1 where the network is defined in a d- 
dimensional Euclidean plane. It was found that the scale 
free behaviour persists above a certain critical value of 
a which depends on the spatial dimensionality Below 
this value of a, the stretched exponential behaviour of 
the degree distribution was again observed. 

We considered a one-dimensional space with periodic 
boundary condition where the nodes occupy the position 
x with < x < 1. Initially we have Too nodes connected 
to each other. First we investigate the scale-free prop- 
erties of the model by studying the degree distribution 
P(k). We vary both a and (3 and observe the behaviour of 
P{k) to obtain a phase diagram. Results for too = to = 1 
and Too = m = 3 showed that the critical behaviour is 
independent of the value of to as in the BA model. We 
noted several interesting features: 

1. In the a — [3 plane, there exists a phase boundary 
along which the network is scale-free. Above this bound- 
ary it shows a gelation-like behaviour as in [16]. Below 
this boundary the degree distribution is stretched expo- 
nential as was observed in [13], [14] and [16]. 

2. Scale free behaviour is observed to occur at the crit- 
ical value (3 C = 1 for all values of a > —0.5. For lower 
values of a it occurs at higher values of (3. For values 
of a < —2.0 the phase boundary is linear given by the 
equation a c + (3 C = 0. 

3. Although the scale- free property is observed along 
the entire phase boundary, there is a difference in the be- 
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haviour of the degree distribution P(k). While P(k) ~ 
fc~ 7 everywhere, 7 ~ 2.7 for a < —0.5 and 7 = 3.0 (as in 
the BA model) for a > -0.5. 

The phase diagram is shown in Fig. 1. We would like to 
emphasise two points from the above observations. First, 
even though the case (3^1 has been studied earlier [16], 
the only point at which scale-free behaviour was observed 
was at f3 c = 1 while here one can get scale-free behaviour 
even at (3 C > 1 by tuning the distance dependence factor. 
Secondly, the exponent 7 = 2.7 ± 0.1 for a < —0.5 may 
not seem to be significantly different numerically from 
the BA value 7 = 3.0 to claim that it belongs to a differ- 
ent universality class. However, as we will discuss later, 
the behaviour of the clustering coefficients are also sig- 
nificantly different here, which will support this claim. 
All the above results were obtained for a network with 
N = 20000 and using 100 different realisations of the 
network. 

We calculated the average shortest path lengths and 
diameter of the model at the phase boundary and found 
that the these two indeed scale logarithmically with the 
number of nodes in the network at the phase boundary 
indicating that the scale-free network also has a small 
diameter. 

The clustering properties of this model are studied in 
detail in an attempt to compare the results with that of 
the real networks. In order to study clustering we kept 
to = mo = 3. Defining the exponents a and b in the 
following way 



and 



C(N) oc N- a 



C(k) oc k~ 



(2) 



(3) 



we find that a and b depend on the values of a and [3. Fig. 
2 shows the behaviour of the clustering coefficients C(N) 
on the critical curve of the phase diagram as a function 
of the number of nodes. We find that for a > —0.5, the 
data is consistent with the behaviour C(N) oc N~ 0J5 . 
The slope of the curves decrease as we go away from 
a = —0.5 to higher negative values indicating that a in- 
creases. This is consistent with the idea that as a is made 
more negative, the nodes get connected to the nearer 
ones making the clustering tendency higher. A curious 
feature of C(N) is that it actually increases with N for 
large negative values of a, e.g., at the critical point cor- 
responding to a = —4.0, a becomes negative. However, 
as the maximum value of C(N) can be unity, we believe 
that a negative value of a indicated that C(N) converges 
to a finite value for (M) — > 00 for large values of a on the 
negative side. 
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FIG. 2. The clustering coefficients as a function of N, the 
number of nodes for different values of a on the phase bound- 
ary (a = —3.0,-2.0,-1.0,0 and 1.0 from top to bottom). 
The gradient in the log-log plot gives the value of a. 

Although the scaling behaviour of C(N) remains same 
for all a > —0.5, calculation of C(N) for a fixed N shows 
that on increasing a the clustering decreases, a result one 
can intuitively guess as for large positive a, the nodes get 
connected to nodes at large distances making the cluster- 
ing tendency lesser. 

Fig. 3 shows the variation of C(k) against k on the 
phase boundary. C(k) is more or less a constant for 
a > —0.5, but for larger negative values of a shows a 
decrease with k. The behaviour of C(k) shows a clear 
power law decay for very large negative values of a. This 
is a feature found in most real- world networks. 
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FIG. 3. The clustering coefficients as a function of k, the 
number of nodes for different values of a on the phase bound- 
ary (a = —3.0,1,5,-1.0 and from top to bottom). The 
gradient in the log-log plot gives the value of b. 

We plot the values of a and b in Fig. 4 at the critical 
points (a c ,f3 c ) as a function of a as we are more inter- 
ested in the role of the spatial distance dependence of the 
network. For a = and (3 = 1, we get the known values 
a = 0.75 and b = 0. For all values of a > —0.5, the 
values of a and b remain the same on the critical phase 
boundary {(3 C — 1) and are the equal to that of the BA 
model. For a < —0.5, the values of a and b are differ- 
ent at different points of the phase boundary. In fact, 
the value of a decreases while b increases towards 1 as a 
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approaches higher negative values. 

We have also studied the behaviour of the clustering 
cocffcients in the regions of the phase diagram where it is 
not scale-free. In the region where there is gelation, C(k) 
shows a power law behaviour again. This is expected 
as most of the nodes get attached to a single node and 
the clustering coefficient decreases as a result. In the re- 
gion where stretched exponential behaviour is observed, 
the clustering coefficient does not show reasonable de- 
pendence on k at large values of k. 
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FIG. 4. The values of a and b on the phase boundary 
(a = a c , (3 — f} c ) as a function of a. 
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To summarise, we have studied a growing network in 
the Euclidean space where the link attachment probabil- 
ity is controlled jointly by two competing factors i.e., the 
preferential attachment and the magnitude of the link 
length. These two factors are tuned by the parameters a 
and (3 as defined in Eqn. (1). A critical boundary in the 
a — (3 phase plane separates the network from its "gel" 
phase to the "stretched exponential" phase. However 
on the boundary between the two phases the network is 
scale-free. Numerical simulations on a one dimensional 
system indicates that on the critical boundary the net- 
work crosses over from a BA universality class (a > —0.5) 
to a new universal scale free behaviour {{a < —0.5). The 
calculation of the exponents a and b for the clustering 
cocffcients defined in Eqns (2) and (3) show that their 
values are non-universal in the region a < —0.5 on the 
phase boundary, with an indication that a converges to 
zero and b converges to unity as a approaches large neg- 
ative values. Thus the network can be tuned to have 
different clustering properties on the phase boundary. 
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In the present model, C(k) oc k~ b with a non-zero value 
of b for a < —0.5 (with (3 — j3 c ) for which the network 
is scale-free and also has a small diameter. The power 
law behaviour of C(k) is obtained as a natural conse- 
quence of (1) without adding further steps in the growth 
process as in the other models considered in recent liter- 
ature. Surprisingly, both the present model and some of 
the other models considered earlier [4,6,7,10] give scale- 
free behaviour as well as C(k) oc k~ b (with 6^0) al- 
though they differ by an important factor - the spatial 
dependence or geographical organisation. Hence, it is 
not possible to guess whether there is any such organi- 
sation present in the network simply by knowing b. The 
real networks with geographical organisation in fact show 
that b = 0, a result we can obtain from the present model 
when the spatial dependence given by a becomes irrele- 
vant and it becomes equivalent to the BA model. Hence 
we conclude that geographical organisation is not the key 
factor responsible for the result 6 = 0. And the result 
b 7^ can be achieved even after incorporating distance 
dependent factors. 

Our present results are for a one dimensional network. 
But as observed in [13], when a^O and (3=1, the two 
dimensional network gives results which are qualitatively 
similar to those obtained in one dimension, we believe 
that in higher dimensions also one would get similar re- 
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